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ABSTRACT 

We describe computational methods for the theoretical study of explicit correlations 
beyond the mean field in excitons confined in semiconductor quantum dots in terms of 
the Auxiliary-Field Monte Carlo (AFMC) method [1]. Using AFMC, the many-body 
problem is formulated as a Feynman path integral at finite temperatures and evaluated 
to numerical precision. This approach is ideally suited for implementation on high- 
performance parallel computers. Our strategy is to generate a set of mean-field states 
via the Hartree-Fock method for use as a basis for the AFMC calculations. We present 
preliminary results. 

INTRODUCTION 

Revolutionary nanofabrication techniques, such as atomic-layer epitaxy, advanced 
lithography, and electrochemical and molecular self-assembly methods, have opened the 
way to fabricating quantum dots, arrays of quantum dots, and other structures that may 
provide the basic building blocks for future nano-scale electro-optic devices [2]. Recent 
experimental studies of optical emission from semiconductor quantum dots have demon- 
strated the importance of many-body effects, such as correlation and exchange, in exci- 
ton dynamics [3]. Intensive basic research is ongoing to accurately simulate the complex 
properties of current and future nanostructure devices using high-performance computers 
(e.g., see [4]). While significant progress has been made in the simulation of the relevant 
independent-particle phenomena, much important work remains in understanding and 
simulating the collective and correlated many-body phenomena observed in nanoscale 
devices at finite temperature [5]. Theoretical research has focused largely on mean-field 
or one-body aspects of the phenomena at zero temperature. Many-body effects, such as 
exciton-exciton interactions, are still lacking a complete theoretical understanding and a 
numerical implementation. 

The focus of this paper is the theoretical description and exact, fully correlated numer- 
ical simulation of many-electron correlations in two-dimensional semiconductor quantum 
dots at finite temperature and exposed to strong magnetic fields. We will approach this 
many-body problem via the Auxiliary-Field Monte-Carlo (AFMC) method [6], which has 
been known since the 1950s [7], but without wide application to fermion systems because 
of difficulties in numerical implementations known as the "fermion-sign" problem. Nu- 
clear physicists have pioneered the modern development of AFMC for the shell-model 



problem [1], obtaining practical solutions to the associated numerical issues. With this 
important work as inspiration, AFMC has been applied in condensed-matter systems 
[8] and atomic and molecular physics [9]. In those cases for which practical solutions 
to the sign problem have been found, the results have been quite revolutionary for the 
theoretical description of the phenomena. No one, to our knowledge, has tried to apply 
the AFMC method to quantum nanostructures. This new effort will provide a numer- 
ically exact treatment of the quantum many-body problem in nanostructures at finite 
temperature. 

In describing the ground-state and low-lying (intraband) excitations of the N-electron 
semiconductor nanostructures, it is often sufficient to restrict consideration to the conduc- 
tion band using the effective- mass approximation [2]. Denoting the operators of creation 
and annihilation of an electron in a single-particle state \a) by c' a and c a , respectively, 
the Hamiltonian of the N-electron system can be written in the occupation-number rep- 
resentation in the form 

Hee = E + W (a(3\V e e\lS)cl&lcsC y , (1) 

« z a/375 

with single-particle energies E e a and matrix elements of the Coulomb interaction between 
electrons (a/3\ V ee \j5). In the following, we will consider nanostructures exposed to strong 
magnetic fields, including the Zeeman spin splitting. For present purposes, we restrict 
ourselves to harmonic-confining potentials. 



AFMC METHOD 



To calculate expectation values in AFMC, we make use of the Euclidian-time many- 
body propagator U = exp(—f3H ee ), where (5 = T _1 is interpreted as the inverse of the 
temperature. For example, we calculate the expectation value of some observable Cl in a 
canonical ensemble, with fixed particle number N, as 

^ = Tr[P N Ub] = Tr N [Un] 
Tr[P N U] Tr N [U] ' 

where the many-body trace is defined as TrX = and the sum is over all many- 

body states of the system. If N is the number operator and Pn = 5(N — N) is the 
projector onto states with A" electrons, the canonical ensemble is defined by Tr N X = 
Ei(i\PNX\i). 

Beyond such static properties, this approach allows one to obtain information about 
the dynamical response of the system. For operators iV and Cl, the response function 
Rq(t) in the canonical ensemble is defined as 

Rn(r) = ^ = <nt(r)n(0)>* , (3) 

where £$(t) = exp(rif ee )iV exp(— rH ee ) is the Euclidian-time Heisenberg operator. In- 
serting complete sets of A^-body eigenstates of H ee |/|)}) with energies E iy f shows 
that 



^(r) = i^e-^|(/|Q|Ore-^/-^, (4) 
*/ 

where Z = J2i ex P( — AE*) is the partition function. Thus, Rq(t) is the Laplace transform 
of the strength function Sq(E): 

/oo 
e- TE S n (E)dE ; (5) 
-oo 

Sn(E) = i£ e-^ | (/|n|<> | 2 5(£ - ^ + . (6) 

fi 

It is important to note that we cannot usually obtain detailed spectroscopic information 
(i.e., energies and wave functions) from AFMC. Rather, we can calculate expectation 
values of operators in the thermodynamic ensembles or the ground state [1]. 

In the AFMC method, two-body interactions in H ee are linearized through the 
Hubbard-Stratonovich (HS) transformation [7]. The difficult many-body evolution U 
is replaced by a superposition of an infinity of tractable one-body evolutions, each in a 
different fluctuating external field, a . Integration over the external fields thus reduces 
the many-body problem to quadrature, which is evaluated stochastically. 

The many-body Hamiltonian can be written schematically as 

H = eO + \v66 , (7) 

where O is a density operator of the form a) a, V is the strength of the two-body interac- 
tion, and e a single-particle energy. In the full problem, there are many such quantities 
with various orbital indices that are summed over, but we omit them here for clarity. 

All of the difficulty arises from the two-body interaction. If H were solely linear in O, 
we would have a one-body quantum system, which is readily dealt with. To linearize the 
evolution, we employ the Gaussian identity 
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dae-iW^e-^ ; h = eO + sVaO . (8) 



Here, h is a one-body operator associated with a c-number field a, and the many-body 
evolution is obtained by integrating the one-body evolution U a = e _/3h over all a with a 
Gaussian weight. The phase s is 1 if V < or % if V > 0. Equation (8) is easily verified 
by completing the square in the exponent of the integrand and then doing the integral. 

With an expression of the form (8), it is straightforward to write observables as the 
ratio of two integrals. For example, the canonical expectation value (2) becomes 

{ )N ~ Jdae-I^Tr N U a ~ I daW a ' W 

where W a = G a TiN U a , G a = e~^ v ^ 2 , and = (Trjv C/ CT fi)/(Trjv Ua)- Thus, the many- 
body observable is the weighted average (weight W a ) of the observable Q a calculated in 
a canonical ensemble involving only the one-body evolution U a . 

An expression of the form (9) has a number of attractive features. First, the problem 
has been reduced to quadrature. Second, all of the quantum mechanics is of the one- 
body variety, which scales simply with the square of the number of single-particle states 



included in the calculation. The price to pay is treating the one-body problem for all 
possible a. 

Since H ee contains many two-body terms that do not commute, the evolution must be 
discretized, i.e., (5 = N t A/3, before applying the HS transformation, i.e., 



Tr N e-^ -> Tr^ [e~ A ^] -> / V[a]G(a)Tr N J] e Af)h ^ , (10) 

J n=l 



where a n are the auxiliary fields at a given imaginary time-step A/5, T>[a] is the measure 
of the integrand, G(a) is a Gaussian in a, and h is a one-body Hamiltonian. Dimensions 
of the integral can reach up to 10 5 for systems of interest, and it is thus natural to use 
Metropolis random walk methods to sample the space. Because the numerical effort 
for AFMC scales only as a low power of the problem size, very large (and hence more 
realistic) calculations are possible. 

The most significant challenge to applying the AFMC method resides in overcoming the 
fermion-sign problem. Algorithmic solutions to these problems are typically not robust 
and are highly dependent on the fermion system to which they are applied. For example, 
systems whose interactions are purely attractive are free of the sign problem [1]. 

For our system of many electrons confined in a quantum dot, the Coulomb forces are 
purely repulsive. In computing expectation values of observables as in Eq. (2), we will 
employ the finite-temperature methods discussed above. At high temperatures, these 
methods will not suffer from a Monte Carlo sign problem, but at lower temperatures, we 
will encounter the sign problem. Its most frequent manifestation occurs when the weight 
functions used for a Metropolis Monte Carlo evaluation of the integrals introduced by 
the Hubbard- Stratonovich transformation lose their positive-definite character [1]. We 
are exploring methods to stabilize the AFMC against the sign problem when the effective 
interactions include terms that are both attractive and repulsive [9,10]. 

Previously, we developed the AFMC method for the nuclear shell-model problem and 
performed state-of-the-art calculations for a wide range of nuclear properties, including 
response functions at finite temperature [1,11,12]. The exciton density in the correlated 
nuclear system has been described by an application of AFMC techniques to the nuclei 
40 Ca and 42 Ca in a very large model space [13]. 

APPLICATION and RESULTS 

We consider the problem of N electrons of effective mass m* in a plane, (x, y), confined 
by an external parabolic potential, V(r) = \m*UQr 2 , and subject to a strong magnetic 

— * 

field B = B e z . We consider the Zeeman splitting but neglect the spin-orbit interaction. 
The Hamiltonian for such a system is 



+ E^^> (ii) 



where the vector potential is A{fi) = (B /2)(—yi,Xi,0), g* = 0.54, e = 12.9, m* = 
0.067m e , and hu = 3meV. 

AFMC calculations require as basic inputs single-particle energies, two-body matrix 
elements of the interaction, and the initial many-body wavefunctions which we obtain 



FIGURE 1. (a) Maximum-density droplet; (b) broken-symmetry phase, "Wigner molecule"; and (c) 
"Wigner crystal" . 

from Hartree-Fock (HF) calculations. The HF method provides physically meaningful 
single-particle states including the exchange interactions exactly [14]. 

Much experimental effort has focused on mapping the magnetic-field dependence of 
the structure of semiconductor quantum dots by measuring the chemical potential via 
capacitance spectroscopy [15]. Cusps and steps in the chemical potential were found to 
clearly separate different ranges of magnetic fields [15,16]. These features were identified 
with phase transitions in the charge density of the quantum dot. For increasing magnetic 
field, all electrons will become spin-polarized initiating the maximum density droplet 
(MDD) phase [15] [see Fig. (la)]. In this phase, the density is constant and homogeneous 
at the maximum value that can be reached in the lowest Landau level. The stability of 
the MDD is determined by a competition between the kinetic and external confinement 
and the Coulomb repulsion between electrons. For increasing magnetic field, the charge- 
density distribution of the droplet reconstructs [18] with a ring of electrons breaking 
off from the MDD phase (see Fig. (lb)). This edge reconstruction has been shown via 
mean-field calculations to result from a rotational symmetry-breaking phase transition 
to a Wigner molecule or Wigner crystal phase [14,17] [see Figs, [lb) and (lc)]. These 
calculations are in good qualitative agreement with recent experimental results [16]. 

While we are in the initial stages of testing and calibrating our AFMC code, we do 
have several interesting initial results. We begin with a canonical ensemble of iV = 6 
electrons in a B = 4 T magnetic field. The corresponding Hartree-Fock solution is shown 
in Fig. (la). We first demonstrate in Figure (2a) the behavior of the Monte Carlo sign ($) 
as a function of inverse temperature j3. Note that beyond approximately /? = 80 eV -1 , 
the sign begins to drop quickly. In order to obtain good evaluations of expectation values 
beyond f3 = 100 eV -1 , we need to calculate many more Monte Carlo samples than here, 
where we have 2560 samples for each (3. Figure (2b) shows the occupation of the various 
Hartree-Fock states at different temperatures. Thermal excitation causes a degradation 
of the spin alignment, and thus an occupation of both spin- up and spin-down levels. 
An interesting near-term calculation will be to find the transition temperature when the 
return of spin polarization occurs. 

We are at the beginning of an exciting and challenging study of the optical properties 
of quantum dots. At the present, we have a partially developed AFMC code with which 
to begin our work. We shall include computations of El dipole-transition operators in 
the code, and will investigate various techniques for overcoming the Monte Carlo sign 
problem for low-temperature studies. 
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FIGURE 2. Left: a) sign problem as a function of inverse temperature; right: b) occupation of the 
HF states, irf with decreasing temperature. 
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